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Abstract 

' Wc present Subaru observations of the newly discovered luminous quasar AKARI J1757+5907, which 

, shows an absorption outflow in its spectrum. The absorption consists of 9 distinct troughs, and our analysis 

focuses on the troughs at ^ —1000 km s~^ for which we can measure accurate column densities of He I*, 
^ Fell and Mgll. We use photoionization models to constrain the ionization parameter, total hydrogen 

I ' column density, and the number density of the outflowing gas. These constraints yield lower limits for the 

2 I distance, mass flow rate and kinetic luminosity for the outflow of 3.7 kpc, 70 Mq yr~^, and 2.0xl0'*'^ ergs 

> . s~^, respectively. Such mass flow rate value can contribute significantly to the metal enrichment of the 

' intra-cluster medium. We find that this moderate velocity outflow is similar to those recently discovered in 

massive post-starburst galaxies. Finally, we describe the scientific potential of future observations targeting 
this object. 

Key words: galaxies: active — quasars: absorption lines — quasars: emission lines — quasars: individual 
. (AKARI-IRC-Vl J1757000+590759) 

,, 1. Introduction provements in analysis methods and in target selection 

(see discussions in Arav et al. 2008; Dunn et al. 2010), we 

Roughly 20% of all quasars exhibit Broad Absorption published several more accurate determinations of these 

Lines (BAL) in their spectrum (Knigge et al. 2008; Hewett quantities (Moe et al. 2009; Dunn et al. 2010; Bautista 

& Foltz 2003), which are indicative of outflows with ve- et al 2010; Arav et al 2010). Here we present a similar 

locities ~ 10'^ — 10^ km s~^. Kinetic energy and mass analysis of an outflow in a luminous overlooked quasar. 
_ ^ emanating from quasars have become key elements in the- Wc use the term "BAL outflow" to designate intrin- 

, oretical modeling of the evolution of supermassive black sic absorption detected in the spectrum of a quasar (i.e., 

' holes (SMBH) and their host galaxies (e.g., Di Matteo et originating from outflowing material in the vicinity of the 

al. 2005; Hopkins et al. 2008), the suppression of cool- AGN, see Barlow & Sargent 1997). The original BAL 

ing flows in clusters (e.g., Ciotti et al. 2010; Briiggen & definition (Weymann et al. 1991) was created to differen- 

Scannapieco 2009), and enrichment of the intra-cluster tiate, in low- resolution spectra, AGN outflow absorption 

and inter-galactic media with metals (e.g., Moll et al. systems from intervening absorbers that do not have a 

2007). Collectively, harnessing quasars' mechanical en- dynamical connection to the AGN and is now physically 

ergy to help in driving the above processes is known as obsolete. Significant number of narrower absorption lines 

"AGN feedback" . show intrinsic natures, i.e., time variability and partial 

To assess the contribution of BAL outflows to AGN coverage of the background light source. The frequency 

feedback scenarios, it is (at least) necessary to determine of intrinsic absorption lines are discussed and summarized 

their mass flow rate {Mout) and kinetic luminosity (Ek). in Ganguly k Brotherton (2008). We therefore use "BAL 

Early attempts to do so were done by de Kool et al. outflow" to designate the physical nature, rather than the 

(2001), de Kool et al. (2002), Hamann et al. (2001), and observational definition, of the phenomenon. 
Wampler, Chugai, Petitjcan (1995). Recently, using im- AKARLIRC-Vl J17570004- 590759 (hereafter 
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AKARI J17574-5907) was discovered during the follow- 

Based on data collected at Subaru Telescope, which is operated observations of AKARI mid-infrarcd (MIR) All-Sky 

by the National Astronomical Observatory of Japan. „ mi • r i . n'. a a t^t r ^ 

Survey. ihe infrared satellite AKARI performed an 
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all-sky survey at 9 and 18 as well as at four far- 
infrared bands (Murakami et al. 2007; Ishihara et al. 
2010). The initial identification of the AKARI MIR 
All-Sky Survey sources involved association with the 
Two Micron All Sky Survey (2MASS) catalog (Skrutskie 
et al. 2006). This search identified some AKARI MIR 
sources with F{9fim)/ F{Ks) > 2 in the high galactic 
latitude (|6| > 20°) after excluding the sample in the 
Large and SmaU Magellanic Clouds. AKARI J1757-h5907 
has a large ratio of mid-IR to near-IR flux density 
{F{9nm)/F{Ks) = 11.1). This mid-IR source is also 
coincident with a bright near-UV and optical source. We 
present photometries of AKARI J1757+5907 in table 1. 
The photographic magnitudes are from USNO-Bl.O cat- 
alog (Monet et al. 2003), and they are converted to g,r, 
and i magnitude using equation 2 of Monet et al. (2003). 
We found the radio source NVSS J175659.82-K590801.5 
(6.5 ± 0.5 mJy at 1.4 GHz) (Condon et al. 1998) is 
coincident with AKARI J17574-5907. The ratio of radio 
(5 GHz) to optical (4400 A) fiux density is 1.4 assuming 
radio spectral index is -0.5. The ratio indicates the 
quasar is radio quiet (Kellermann et al. 1989). 

The follow-up spectroscopy using KPNO 2.1m telescope 
revealed that AKARI J1757-I-5907 is a z = 0.615 quasar 
that shows He I* and Mg II absorption lines as well as H/3 
and strong Fc II emission lines (Toba et al, in prepara- 
tion). The spectrum resembles the one of QSO 2359- 
1241 (Brotherton et al. 2001; Arav et al. 2001). Both 
quasars show rare He I* absorption as well as Mgll ab- 
sorption. They also have redder continua and strong 
Fell emission lines. The high resolution spectroscopy of 
QSO 2359-1241 by Arav et al. (2001) revealed Fell ab- 
sorption lines, Thus Fe II absorption lines are expected 
to exist in AKARI J17574-5907. Following Korista et al. 
(2008), by measuring the ionic column densities {Nion) of 
Fell and He I*, we can constrain the ionization parame- 
ter {Uh) and hydrogen column density (Nh)- The high 
brightness of this quasar permits us to do high resolution 
spectroscopy of He I* absorption lines and search for Fe II 
resonance and excited state absorption lines. The Nion 
ratio of Fell* to FeIl(E=0) yields the hydrogen number 
density, which in turn yields the distance of outflowing 
gas from the central source, Ek and Mout by using the 
number density, Uh and Nh- 

The plan of this paper is as follows. In § 2 we describe 
the observations and data reduction. In § 3 we deter- 
mine the redshift of the object. The outflow absorption 
troughs are discussed in § 4 and the spectral energy dis- 
tribution in § 5. In § 6 we describe our photoionization 
modeling, and in § 7 the resultant determination of mass 
flow rate and kinetic luminosity for the outflow. In § 8 
we discuss our results and in § 9 we describe the scien- 
tific potential of additional Subaru/HDS and HST/COS 
observations targeting this object. Throughout this pa- 
per we assume Hq = 70 km s~^ Mpc~^, il„i = 0.3, and 
= 0.7. Note that wavelengths of any transition in this 
paper are ones in vacuum, and the observed wavelength 
scale is converted to one in vacuum. 



2. Observations and Data Reduction 

2.1. high resolution spectroscopy 

The high resolution spectroscopy of 
AKARI J1757-t-5907 (R.A.=17:57:00.24, 

Decl.=-t-59:08:00.3 (J2000.0)) was done with the High 
Dispersion Spectrograph (HDS; Noguchi et al. 2002) 
attached to the Subaru 8.2 m telescope (lye et al. 2004) 
on 2010 June 17 (UT). The weather conditions were poor 
with thick cirrus and the seeing was unstable (>1"0). 
The slit width was set to be 1"0. The HDS setting 
was Yd covering between 4054 and 6696 A. This results 
in a resolving power oi R ^ 36000. The binning was 
2 (spatial direction) x 4 (dispersion direction). We 
obtained eight exposures of 1800 s each, however, one of 
them is unusable due to low signal-to-noise ratio. 

The data were reduced using IRAF^ for the standard 
procedures of overscan subtraction, dark subtraction, cos- 
mic ray removal and flat-fielding, where wavelength cal- 
ibration was performed using the Th-Ar lamp, the rms 
wavelength calibration error is 0.011 - 0.013 A. The one- 
dimensional spectra were extracted from each exposure. 
Heliocentric correction was applied. After that, all the 
spectral exposures were combined, and all orders were 
connected to one spectrum. We normalized the spectrum 
using spline fits. Finally we converted to vacuum wave- 
length scale. 

2.2. Spectrophotometry 

The low resolution spectrophotometry of 
AKARI J1757+5907 was done on 2010 June 30 (UT) with 
FOCAS (Kashikawa et al. 2002) attached to the Subaru 
telescope. The new fuUy-depleted-type CCDs developed 
by NAOJ/ATC and fabricated by Hamamatsu Photonics 
K. K. were installed and commissioned at that time. We 
obtained six spectra of 5 minutes integration under a clear 
condition and good seeing (0."6 - 0"9 ). The 2."0 width 
slit was used for spectrophotometry purpose. We used 
two configurations: the R300 grism with the 058 filter 
('red') and the B300 grism without any order-cut filters 
('blue'). The first three 300s spectra were the 'red' which 
covers between 5700 A and 10200 A. The last three 300s 
spectra were the 'blue' which have an uncontaminated 
range between 3500 A and 7000 A. The atmospheric 
dispersion corrector was used. The slit position angle 
was 0°, and the binning was 2 (spatial direction) x 1 
(dispersion direction) The spectrophotomctric standard 
star BD+28°4211 was observed for sensitivity calibration. 

The data were reduced using IRAF for the stan- 
dard procedures of bias subtraction, wavelength cali- 
bration, and sky subtraction, except for fiat-fielding. 
AKARI J1757-I-5907 is so bright that its counts on 
the CCD are comparable to the flat frames, and are 
much higher at shorter wavelength (< 4000 A). The flat- 
fielding procedure significantly reduced its signal-to-noise 

^ IRAF is distributed by the National Optical Astronomy 
Observatory, wliicli is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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ratio. We therefore skipped the flat-fielding procedure. 
Wavelength calibration was performed using OH night sky 
emission lines for the red spectra and the Cu-Ar lamp for 
the blue ones. The rms error of wavelength calibration is 
0.2 A. The seeing was much smaller than the slit width, 
thus the resolution was determined by the seeing disk size. 
The resulting He I* 3889 absorption line is 13 A width at 
6260 A. This value corresponds to the resolving power of 
480 (~ 620 km s~^), which is similar to the resolution ob- 
tained with the 0."8 width slit, and is consistent with the 
seeing size during our observations. The sensitivity cali- 
bration was performed as a function of wavelength. The 
flux of the blue and red spectra at the same wavelength 
agree within 1.6%. The foreground Galactic extinction 
of E{B -V)= 0.043 mag (Schlegel, Finkbeiner, & Davis 
1998) was corrected. 

3. Results of spectrophotometry 

Figure 1 displays the low-resolution optical spectrum 
of AKARI J1757-I-5907. The spectrum shows strong ab- 
sorption lines of Mg II and He I* as well as emission lines of 
Mgll, Fell, H7, H/3, [O III]. In order to determine the sys- 
temic redshift of the quasar, we deblend the [O III] emis- 
sion lines from the H/3 and Fell emission lines. As seen 
in flgurc 1, AKARI J17574-5907 has strong Fell emission 
lines. However, the Fell emission template around H/3 
derived from the spectrum of I Zw 1 (Aoki, Kawaguchi, 
& Ohta 2005) is not a good match to the Fell emission 
in AKARI J1757-f5907 (figure 2a). The intensity ratios 
among Fe II emission lines are different between these ob- 
jects. Thus, we cannot use the Fell emission template 
derived form the spectrum of I Zw 1. 

Instead we fit the spectrum between 7662 A and 8204 A 
with a combination of Fell AA4925,5020, H/3, and [Olll] 
after subtraction of a power-law continuum. This power- 
law continuum is constructed by fitting at 6526-6544 A 
and 8842-8891 A. The H/3 emission line is fitted with a 
combination of three Gaussians. We fit the [Olll] dou- 
blet AA4960. 3, 5008.2 with two sets of two Gaussians. The 
width and redshift arc assumed to be the same for each set, 
and the intensity ratio is fixed to be 3.0. Fell AA4925,5020 
are modeled by two Lorentzians with the same width and 
redshift. Their redshift is fixed to be the same as the 
strongest Gaussian component of the H/3. The result of 
fitting is shown in figure 2b. We need a separate "red 
wing" of H/3 to get a satisfactory fit. This component 
may be [O III] emission line. The derived redshifts and 
FWHMs are tabulated in table 2. The FWHM is cor- 
rected for the instrumental broadening by using the simple 
assumption: FWHMtruc = (FWHM^i^^ - FWHMf^.Ji/a^ 
where FWHMobs is the observed FWHM of the line and 
FWHMinst is an instrumental FWHM (620 km s'^). The 
redshift of the red component of [Olll] is 0.6150 ± 0.0001. 
The rest equivalent width of [O III] including both compo- 
nents is 9.9±0.3 A. This value is consistent with the [O III] 
strength in majority (> 50 %) of Mgll BAL quasars re- 
ported by Zhang et al. (2010). We also detect weak [Oil] 
emission at 6023.4 A, which corresponds to a redshift of 



0.6155. We thus adopt 0.61525 ± 0.00025 as the systemic 
redshift of the quasar. The blue component of [O III] is 
blueshifted by 980 km s~-^ relative to the systemic red- 
shift. 

4. Absorption lines 

In the HDS spectrum of AKARI J1757-h5907, the ab- 
sorption lines are not heavily blended. The identifica- 
tion is thus a straightforward task. We identify Mgll 
AA2976,2803, He I* A2945, A3188, A3889, Fell A2600, 
A2586 as well as weaker absorption troughs from Mgl 
A2852, He I* A2829, and Call AA3934,3969. The strong 
absorption lines such as the Mgll doublet and HeI*A3889 
clearly show 9 distinct troughs, which span a velocity 
range from -660 to -1520 km s"^ The trough at -1000 
km s~^ has the same outfiow velocity as the blue com- 
ponent of the [Olll] emission line. We show the ab- 
sorption troughs from Mgll AA2976,2803, He I* A2945, 
A3188, A3889, and Fell A2600, A2586 in figure 3. We 
do not detect an Fell A2612 absorption trough from the 
E=385 cm~^ excited level, which has the largest oscillator 
strength of the Fell* lines from this energy level, in our 
spectral range. We show the spectral region of Fell A2612 
in the lower panel of figure 3. 

4-.1. Column Density Determinations 

Both the Mgll and He I* troughs span the full veloc- 
ity range from —660 to —1520 km s~^ and appear in all 
9 distinct troughs (see Figure 3). There is self blend- 
ing in the Mgll troughs at the velocity extremes, which 
does not affect the majority of the troughs. The He I* 
and Fe II troughs are free of any self blending. Of the 9 
troughs, only three have corresponding absorption in Fe II. 
These are in the range of —1050 to —800 km s~^ (see 
Figure 3). Thus, we concentrate on this velocity range for 
our measurements as the best photoionization constraints 
are achieved by contrasting He I* and Fe II column densi- 
ties (see section 6; Korista et al. 2008; Arav et al. 2010). 

To determine the ionic column densities, we rebin the 
data to 10 km and use the velocity dependent ap- 
parent optical depth (AOD), covering factor (C(v)), and 
power-law fitting methods from Dunn et al. (2010; meth- 
ods 1, 2, and 3 in their subsection 3.2) across the range of 
-1050 to -800 km s~\ 

The AOD methods assumes that the emission source is 
completely and homogeneously covered by the absorber, 
so that the optical depth (t(u)) at a given velocity is re- 
lated to the normalized intensity via: I{v) = cxp(— r(u)). 
The covering factor method assumes that at a given ve- 
locity, a fraction C{v) of the emission source is covered 
with a constant value of optical depth, while the rest 
of the source is uncovered. In order to for both t{v) 
and C{v) we use use at least two absorption lines from 
the same energy level of the same ion, and solve for 
I{v)j = 1 — C{v) + C {v) exp{—T{v)j) , where Ij is the nor- 
malized intensity of the absorption due to the j tran- 
sition in the same energy level. The ratio of different 
T{v)j are known from atomic physics and therefor the 
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set of equations is solvable. The power-law model as- 
sumes that the absorption gas inhomogeneously covers 
the background source. The optical depth is described 
by Ty{x) = Tmax{v)x°' , where x is the spatial dimension in 
the plane of the sky, a is the power law distribution index 
and Tmax is the highest value of r at a given velocity. In 
case of the power law model, r is averaged over the spatial 
dimension x. 

In order to convert t in each model to column density 
we use 

N{v) = 3.8 X 10"-^r(i;)(cm"2]^m"^s), 

where A is the wavelength of the line in A, / is the oscilla- 
tor strength. We use the oscillator strengths from Fuhr & 
Wiese (2006) for the Fell lines and values from the NIST 
Atomic Spectra Database (2010) for Mgll and He I*. 

In Figure 4, we show the He I* troughs and column den- 
sity determination from the three lines of He I* present in 
the HDS spectrum. The He I* lines are well separated, 
have no self blending, and we obtain consistent results 
with all three column density extraction methods. There 
are two velocity points where the C{v) becomes nonphys- 
ical (i.e., negative or larger than 1.0), near —930 km s^^ 
and at velocities lower than —850 km s~^. This occurs 
because the two weaker lines are very shallow at these ve- 
locities and are thus consistent with the continuum level 
and dominated by the noise. 

Wc show in Figure 5 the trough profiles and column 
density results for Mgll. Both the red and the blue dou- 
blet troughs are quite deep and therefore their level of 
saturation is model dependent. The C{v) solution sug- 
gests a small level of saturation (only 40% larger column 
density than the AOD estimate), while the power-law re- 
sult is three times higher. This is an inherent feature 
of the absorption models where in order to fit deep dou- 
blet troughs the power-law model requires a much greater 
column density than the C{v) solution (see the case of 
the O VI doublet in the spectrum of Mrk 279, Arav et al. 
2005). Since we do not have data for troughs from ad- 
ditional Mgll lines, we cannot determine which model is 
more physical. 

There is a possibility that a narrow Mg II emission line 
fills the trough. As already noted early in this section, 
the blue component of [Olll] has the velocity of -1000 km 
s~^. The Mgll emission line from the same gas probably 
exists. We estimate the flux of Mgll emission line using 
the flux ratio of [O III] A5008 to Mg II in Seyfert 2 galaxies, 
NGC 1068 (Kraemer et al. 1998) and Mrk 3 (Colhns et 
al. 2005). The ratios vary along the physical position in 
the galaxies between 0.03 and 0.15, and average of extinc- 
tion corrected ratios are 0.09, and 0.07 in NGC 1068 and 
Mrk 3, respectively. We also assume a 1:2 ratio for the 
Mgll doublet, and gaussians of the same width {a = 205 
km s~^) as the [Olll] emission line. This width corre- 
sponds to 3.0 A at the Mgll observed wavelength of 4500 
A. The HDS spectrum before normalization is scaled to 
the low-resolution spectrum. The expected height of the 
Mgll emission line will be ~ 13 % and ~ 4% of the residual 



intensity at Mgll A2796 and A2803 trough, respectively. 
Thus, the real residual intensity may be smaller than the 
observed one. Therefore, we conclude that the Mgll col- 
umn density can be much larger than the AOD or C{v) 
determined values. 

Finally, we show the result for Fell troughs in Figure 6. 
Unlike Mgll and He I*, the Fell troughs are both shallow 
and in a much lower signal-to-noise region of the spec- 
trum (towards the short wavelength end of the detector) . 
Due to this, we find that the weaker Fell A2587 line is 
only detected across three velocity bins. Using both the 
Fell A2600 and the Fell A2587 lines we calculate Npeii 
for both C{v) and power-law methods for the three bins 
and include them in the integrated total. Due to a lack 
of detection of the A2587 line, we use the A2600 line to 
calculate the column density from the AOD method for 
the remaining points. We also checked the column density 
for Fell using the covering factor of Mgll. The column 
density calculated in this fashion changed only by ^10% 
as Mgll nearly fully covers the source (C(w) saO.9) in this 
velocity range. 

4-2. Column Density Limits on Fell* Excited State 
Lines 

In order to help constrain the photoionization models 
in section 6, we estimate the column density limits for the 
Fell* A2612 and A2757 lines from the 385 and 7955 cmT^ 
energy levels, respectively. Neither line shows a detectable 
trough in the data. Therefore, we can use the trough pro- 
file of Fe II A2600 to determine the relative optical depth 
and column density of these two energy levels (see Section 
3.3 of Dunn et al. 2010). We find upper limits on the ionic 
column densities of (3.7 and 0.5) xlO^^ cm~^ for the 385 
and 7955 cm~^ energy levels, respectively. 

5. Determination of the Spectral Energy 
Distribution 

Our FOCAS spectrophotometry data and the GALEX 
photometry clearly show the flux drops at shorter wave- 
lengths (< 2000 A in the rest frame) and indicate redden- 
ing by dust (figure 7). We must consider extinction for 
deriving the intrinsic spectral energy distribution (SED). 
First, we measured the continuum fiux at four points 
where there are less Fell emission contaminates. These 
four points and the GALEX photometry points are then 
shifted to the rest frame. We fit a reddened power-law 
to the continuum points. Wc adopted the index of the 
power-law continuum (a; ccv"') of the LBQS compos- 
ite, -0.36 measured by Vanden Berk et al. (2001). We 
adopt the SMC-type extinction law. The best fit gives 
us the color excess E{B — V) of 0.18. The best fit of the 
reddened power-law continuum is shown in figure 7. 

To estimate the distance to the outflow from the central 
source (i?), we need to determine the flux of hydrogen 
ionizing photons that irradiates the absorber (see equation 
[1] below). Using the Mathews & Ferland (1987) SED, 
reddened to match the observed spectrum, we find that 
number of hydrogen ionizing photons emitted per second 
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by the reddened central source {Qh) is 2.2x10° photons 
s^^. Here we assumed that the reddening occur between 
the central source and the outflow, as is the case where 
the photons are attenuated by the edge of the putative 
AGN obscuring torus (see full discussion in Dunn et al. 
2010). This assumption will also give us smaller values 
for the inferred R and therefore conservative lower limits 
for Alout and Ek- The LboI for the dereddened, intrinsic 
spectrum is 3.7 x 10^^ erg s~^. 

6. Photoionization Modeling 

Through photoionization modeling, reliable measure- 
ments of He I* and Fe II column densities provide accurate 
constraints on the total hydrogen column density, Nh, 
and the hydrogen ionization parameter, 

Qh 



Uh = 



(1) 



where R is the distance from the central source, nn is 
the total hydrogen number density, and c is the speed 
of light. We use version cOS.OO of the spectral synthesis 
code Cloudy, last described by Ferland et al. (1998), to 
model a plane-parallel slab of gas with constant hydrogen 
number density irradiated by a source continuum. We 
focus on the kinematic components spanning a velocity 
range from —800 km s~^ to —1050 km s~^ where we de- 
tect Fell(O). In this velocity range, measurements of up- 
per limits on the FeIl(E=385 cm^^) and the FeIl(E=7955 
cm~^) yield upper limits of electron number density rie, 
lO'^'^ and 10^-^ cm~^, respectively. We adopt the conser- 
vative value of an upper limit of < 10^'^ cm~^, and 
assume tih ~ Tif,, which is valid within the ionized zone 
we are discussing. While the electron number density is 
well constrained from above, there are no diagnostics for 
a lower limit on rie in the data. 

We begin investigation of the parameter space by using 
Cloudy's optimization mode to determine Nh and Uh 
for nH = lO'^'* cm~'^. The parameters Nh and Uh are 
varied and ionic column densities are computed for each 
set of parameters. Best fit values are determined by 
minimization for given tolerances in the measured ionic 
column densities. We adopt the measured ionic column 
densities determined by the partial covering method, and 
optimize to Npe ii and Nne i*, since these are the more 
robust measurements. The measured and model predicted 
column densities are presented in table 3. For solar abun- 
dances and the MF87 SED, as implemented by Cloudy^, 
we find log Uh = -2.15 and log Nh = 20.82 yield good 
fits to the column densities of He I* and Fe II while Nmq n 
is overpredicted. However, as discussed in section 4, the 
Mg II troughs may be more saturated than the partial cov- 
ering model suggests. A hydrogen ionization front, which 
we define as the position at which half of the total hydro- 
gen is neutral (approximated by Nh = IQ'^^Uh), does not 
form in this solution although we are very close to it with 
logiNH/Un) = 22.96. 



Since the data do not provide a lower limit to the elec- 
tron number density, we find other valid solutions by re- 
ducing Uh- When the hydrogen number density is re- 
duced, the He I* population drops. This is due to the 
fact that He I* is populated by recombination of He II and 
the number of recombinations per unit time depends lin- 
early on rie. Therefore, Nh must increase in order to 
provide enough He I* to be consistent with the measured 
value. Fell becomes dominant near the hydrogen ioniza- 
tion front, while He I* drops off drastically at the front. 
Thus, solutions in this region of the slab have Uh fixed 
by Nhc I* and the ratio Nh/Uh fixed by Npe u- Due 
to the tight correlation of Nh and Uh, all valid solutions 
have a similar ratio laying (nearly) on a straight line in 
the Nh-Uh plane for number densities greater than ^ 100 
cm~^. At these densities, the slabs do not form a hydro- 
gen ionization front. As we go to densities lower than 
^ 100 cm~'^, a front forms, behind which Fell and Mgll 
increase linearly with Nh- 

In order to determine what effects choice of SED may 
have on Nh and Uh, we compare results of the MF87 SED 
with results of a softer SED. The soft SED has an opti- 
cal to X-ray spectral index aox = — 1.5 compared to the 
MF87 SED with Uox = -1.4 (with the convention F^^i^") 
and was generated using the Cloudy command agn 375000 
-1.50 -0.125 -1.00, where the numbers are the tempera- 
ture of the UV bump, aox, ctuv, ctx, respectively. We 
find the resulting Nh and Uh are nearly identical and 
conclude that changing the SED only affects the energet- 
ics through Qh, a finding consistent with the analysis of 
QSO 1044-^3656 reported in Arav et al. (2010). 

Another assumption in our models is solar abundances. 
To check the sensitivity of our results to metallicity 
changes, we use the abundances in table 2 of Ballero et 
al. (2008) for metahicity Z=4.23 with the MF87 SED and 
lognn = 3.8. While helium abundances are expected to 
increase with oxygen abundances (e.g., Olive & Scully 
(1996)), the amount of increase varies for different galax- 
ies. Therefore, to be conservative, we increase the helium 
abundance significantly to 15% above solar. We find that 
logUn is approximately 0.02 dex lower and logNn is ap- 
proximately 0.15 dex lower for the increased metallicity 
model with logiNH/Un) = 22.83. We discuss the effect of 
these changes on the energetics of the outflow in the next 
section. 

7. Energetics of the Outflow 

Of particular interest for any outflow are its mass 
(Mout), the average mass flow rate (Mout) and mechani- 
cal work output or kinetic luminosity (i?fc). Assuming the 
outflow is in a form of a thin partial shell moving with 
a constant radial velocity (v), at a distance R from the 
source, the mass of the outflow is: 



■ Att ^uipVtR^ N H , 



(2) 



2 This SED differs from the MF87 SED by the addition of a sub- 
millimeter break at 10 /xm. 



where Nh is the total column density of hydrogen, nip is 
the mass of a proton, ^=1.4 is the plasma's mean molec- 
ular weight per proton, and VL is the fraction of the shell 
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occupied by the outflow. The average mass flow rate is 
given by dividing the outflowing mass by the dynamical 
time scale of the outflow R/v (see full discussion in Arav 
et al 2010), therefore 

Mont ~ = An^impnRNHV (3) 

K/ V 

and Ek = ^Moutv"^ ^ 2-k ^rupVlRN hv^ . (4) 

For the troughs we consider in AKARI J1757+5907, 
the median velocity of the system is —970 km s~^. We 
assume that = 0.2, which is the percentage of quasars 
showing BALs in their spectrum (see discussion in Dunn 
et al. 2010). 

To determine the distance, we solve equation (1) for i?, 
which depends on Uh and nn- The lack of Fe II* detection 
yields an upper limit of rin < 10'^'^ cm~^ (see Section 6), 
and therefore, as shown below, a lower limit on R. 

For the range 10^'^ < nn < 10"^ *, our photoionization 
solutions obey the relationships Nh oc Uh and Uh oc n^", 
where 0.4 < a < 1 (a decreases as nn increases). The first 
relationship arises from the requirement of being close to a 
hydrogen ionization front, and the second is due to the de- 
creasing electron population at the He I* mcta stable level 
for lower number densities. Therefore, equation (1) yields 
R (X ^ ^J^J'Q ^j^e solar abundances photoion- 

ization solution values from Section 6 (logJ7/f = —2.15, 
logNH = 20.82 cm-2 for nn = 10^'^ cm^^) ^nd Qh from 
Section 5. Inserting these values into equation (1), pro- 
vides a lower limit on the distance, R > 3.7 kpc for 
Uh = 10^'^ cm~^, which only increases to 6.6 kpc for 
Uh = 10^'^(cm~'^) due to the weak dependence of R on 
Uh- From equation (3) we observe that Mout and Ek de- 
pend linearly on the product RNh, which is proportional 
to rijj"'^'^"^^'^ . Therefore, the upper limit for nn provides 
lower limit to Mout and Ek ■ Using the Nh and R derived 
for TiH = lO'^'^ cm~'^, we find lower limits of Mout > 70 
^^0.2 yr~^ and Ek > 2.0 xlO^^ O0.2 ergs s~\ where 
SI0.2 = ^^/0.2. 

In the previous section, we discussed changes in Nh and 
Uh due to metallicity and SED changes. Using the Ballero 
et al. (2008) abundances for Z/Zo=4.23 reduces the mass 
flow rate by ^ 30%. Changing to the soft SED described 
in the previous section results in very small changes in Nh 
and Uh, but Qh increases by a factor of ^ 2, increasing 
the mass flow rate by a factor of ^ -\/2- 

The mass of the black hole (Mbh) of 
AKARI J1757+5907 is derived to be 4.0 x IO^Mq 
based on the width of the H/3 emission line of cr = 2160 
km s~^ and the dereddened optical luminosity AL(5100A) 
of 3.8 X 10^^ erg s~^. We use the formula in Bennert et 
al. (2010) based on the calibrations of broad-line region 
size-luminosity relation (Bentz et al. 2006) and the virial 
coefficient taken from Onken et al. (2004). The derived 
log Lhoi/LEdd = —0.13 is used for calculation of mass 
accretion rate (Mace) based on the accretion disk model 
by Kawaguchi (2003), which takes into account the effects 
of electron scattering (opacity and disk Comptonization) 



and the relativistic effects. The Mace is 110 Mq yr 
which is similar to the lower limit of Mout ■ 

Wc note that the ft ~ 0.2 wc use is based the percentage 
of quasars showing C IV BALs. In a recent work, Dai et 
al. (2010) have shown that in the near infra-red surveys 
low ionization BALs (LoBALs) fraction is 4%, consider- 
ably higher than deduced from optical surveys (probably 
due to obscuration effects in the optical band). A more 
appropriate comparison in our case is to include some- 
what narrower outflows with 1000 km s'^ < Av < 2000 
km s~^. For these LoBALs based on "Absorption index" 
(Trump et al. 2006), Dai et al. (2010) find a 7.2% fraction. 
The frequency of He I* outflows is much less known. The 
strongest He I* line in the optical (A3889) is shifted out- 
side the optical range for objects where we can detect C IV 
A1550 from the ground (e.g., z=1.5 for SDSS spectra), so 
a meaningful census of these outflows is difficult to come 
by. Anccdotally, wc flnd that in most cases where we de- 
tect Fe II absorption trough, we also detect He I* troughs 
provided wc have a clear spectral coverage of the latter 
(e.g., Arav et al. 2008; Arav et al. 2010). 

We also point out that AKARI J1757-I-5907 as weU as 
QSO 2359-1241 do not have measurable Fell absorption 
in their low resolution spectra. Also, their outflow veloci- 
ties (^ 1200 km s~^) and widths of troughs are moderate. 
Quasars with similar moderate outflows are more numer- 
ous than extreme FeLoBALs (e.g., SDSS J0318-0600, Hah 
et al. (2002)), and as we show here, can have similar mass 
flow rate and kinetic luminosity as the more extreme ones. 
This fact suggests the mass flow rate and kinetic luminos- 
ity values found here are more common among quasars, 
than judged by the rarity of extreme FeLoBALs. If we 
assume as a conservative limit that the of the outflow 
seen in AKARI J1757+5907is similar to the 7.2% LoBALs 
fraction found by Dai et al. (2010) then fig. 2 = 0.36, which 
will reduce the values for mass flow rate and kinetic lumi- 
nosity accordingly. 

8. Discussion 

In Table 4 we show our Mout and Ek determinations in 
quasar BAL outflows to date (the older, less accurate, ones 
are shown in table 10 of Dunn et al 2010). While the lower 
limit on Ek for the AKARI J1757-I-5907 outflow is rather 
low for AGN feedback purposes, we note that the corre- 
sponding Mout value is large enough to yield a significant 
contribution for the metal enrichment of the intra-cluster 
medium around the parent galaxy (see Hallman and Arav 
2010 [ApJ submitted]) 

We also note that this moderate velocity outflow is 
similar to those discovered in post-starburst galaxies at 
z ~ 0.6 (Tremonti, Moustakas, Diamond- Stanic 2007). 
Those galaxies show an outflow with Mgll absorption 
and a velocity of 1000 km s^^ . Their stellar masses are 
as high as (0.7 — 4.8) x IO^^Mq (calculated from stellar 
population synthesis modeling flt to their spectra). For 
comparison, the Mbh of AKARI J1757-I-5907 is 4.0 x 
IO^Mq, therefore based on Msn-Mbuige relation (Haring 
& Rix 2004), the bulge of its host galaxy is estimated 
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at Mbuige^-^ X lO^^Af0. In addition, our derived Nh is 
similar to the one Tremonti, Moustakas, Diamond-Stanic 
(2007) derived for outflow in the massive post-starburst 
galaxies {Nh — 2x 10^° cm s^^). They also pointed out 
that AGNs probably exist in those post-starburst galaxies 
because they get better fit to the spectra with a feature- 
less blue power-law continuum and they detected high- 
excited emission-lines such as [Olll] and [NeV]. Thus, 
the outflows seen in the massive post-starburst galaxies 
may be same phenomena as outflows in quasars. These 
facts may be one of the indications that the outflow is a 
common phenomenon among massive galaxies. Tremonti, 
Moustakas, Diamond-Stanic (2007) had to make several 
assumptions to derive the physical quantities of the out- 
flow because high resolution spectroscopy for such faint 
targets is difficult and non-detection of the important di- 
agnostic absorption lines from Fe II (both ground and ex- 
cited states) and He I*. In contrast, our high-resolution 
spectroscopy of BAL outflows in quasars, can constrain 
the total hydrogen column density and the distance of 
outflowing gas from the nucleus, which yield less model 
dependent estimates for Mout and Ek- 

9. Scientific gains with additional observations 
9.1. Additional Subaru HDS data 

We expect data with a much higher signal-to-noise ratio 
using Subaru/HDS under good weather conditions. A half 
hour exposure should have a signal-to-noise ratio of 35 
per pixel at 3800 A where Fe II* A2396 is shifted to. Our 
current upper limit to Fc II* A2612 is derived from the data 
with signal-to-noise ratio of 30 per pixel. Fell* A2396 is 
2.3 times stronger absorption than Fell* A2612. With 
two hours (4 x 1800 sec) of exposure time, we expect five 
times more strict upper limit of the column density for 
the excited state of Fell (385 cm~^). Such an upper limit 
would translate to a hydrogen number density of nn < 
10^-^ cm~'^. For this value of tih and our best model at 
that density {logUH ~ -1.66 and logNn = 21.32), we find 
R > 5.9 kpc, Ek > 1-2 xl044 fi^ 2 ergs s'^ and Mout > 
340 r2o.2 Mq yr~^. It is of course possible that we'll be 
able to detect Fe II* A2396 absorption, which will give us a 
determination of and therefore measurements (instead 
of lower limits) for Mout and E^- 

In this paper our analysis focused on the troughs for 
which we can measure Fell column density. The other 
six troughs have only upper limits for Fell column den- 
sity given the current signal-to-noise ratio. In order to 
study the relationship between the troughs, and obtain 
the integrated values of Mout and Ek for the full outflow, 
it is important to measure Fell column density in all the 
troughs. Combined with our other measurements, this 
will constrain the ionization parameter and total hydro- 
gen column density in each individual trough. We detect 
He I* absorption from 3 different lines in almost all the 
troughs (figure 3). The current data at He I* 2945,3188 
have a signal-to-noise-ratio of 65 per pixel. In order to get 
similar or better signal-to-noisc ratio at Fe II resonance 
lines, 2383,2586, and 2600, we will need two additional 



hours of HDS integration. 

9.2. Imaging of outflowing gas 

The blue component of [O III] emission line has the 
same outflow velocity to the trough at ~ 1000 km s^^. 
Furthermore, the derived density and ionization param- 
eter values for the trough are typical for the narrow-line 
regions of AGNs (see, e.g.. Groves, Dopita, Sutherland 
2004). The scale of 3.7 kpc in AKARI J1757-h5907 cor- 
responds to 0"55. The extended nebular gas associated 
with the outflow can be observed by the [O III] emission 
line by HST or the [S III] A9533 emission line by integral 
filed spectroscopy coupled with adaptive optics observa- 
tions from large ground-based telescopes. Currently, the 
solid angle subtended by the outflow is the parameter with 
the largest uncertainty. Imaging the outflowing gas will 
directly determine this parameter. 

9.3. HST Cosmic Origins Spectrograph (COS) observa- 
tions 

The near UV side of COS offers two important 
gains for the scientiflc investigation of the outflow in 
AKARI J1757-f 5907. First, using the G230L grating we 
can cover the full range of 1333-2800 A (observed frame) 
with resolution of ^ 3000 and obtain data with S /N ^ 30 
per resolution element using a total of six HST orbits. 
Such data will allow us to connect the low ionization ab- 
sorption studied in this paper with the higher ionization 
phase seen in Si IV, CiV, and Nv. In addition, this spec- 
tral range covers 3 pairs of Sill/Sill* lines that can deter- 
mine number density considerably lower than is possible 
with the lines from the Fell E=385 cm^^ lines (essen- 
tially the critical density of the Si II* level is an order of 
magnitude lower than that of the Fell E=385). In the 
unlikely event that we will detect Sill, but not Sill* ab- 
sorption, which means that the gas density is substan- 
tially below the critical density for the Si II* level, we will 
have C Il/C II* transitions covered that can determine the 
number density down to ~ 10 cm~'^. These diagnostics 
practically guarantee that we will be able to determine 
He and therefore the distance of the outflow, as long as 
R < 100 kpc. 

In addition we can target the strongest pair of Si Il/Si II* 
fines (1260 A, 1265 A) for a higher resolution (i?- 20,000) 
in order to obtain a fully resolved trough where we can 
obtain function of velocity for all 9 outflow com- 

ponents. With 6 HST orbits using the COS 225M grating 
we can obtain high enough S/N to determine for 25-50 
resolved velocity points across the full width of the out- 
flow for a large dynamical range in rie (100 < rie < 3000 
cm^"^). This will allow for a sensitive tomography of the 
outflow and precise determination of the distance to each 
kinematic component. 
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pecially T.-S. Pyo, and A. Tajitsu, for their assistance 
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taking FOCAS spectra kindly and providing to us. We 
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No. 1 AKARI J1757+5907 



Rest wavelength (A) 

3000 4000 5000 



6 X 
o ^ 




4000 6000 8000 

Observed wavelength (A) 

Fig. 1. Low resolution spectrum of AKARI J1757+5907. 
Ordinate is a flux density corrected for the Milky Way ex- 
tinction in units of erg s~^ cm~^ and abscissa is the 
observed wavelength in angstrom. The rest wavelength is 
given along the top axis. The hatched regions indicate the 
place of telluric absorption. 



Velocity (km s ) 



Fig. 3. Outflow troughs in AKARI J1757-I-5907. Ordinate 
is a normalized flux density, and abscissa is velocity from the 
systemic redshift {z = 0.61525). Blended parts are denoted as 
dotted spectra. The dashed vertical line indicates the position 
of the blue component of [O III] emission line. The velocity of 
that [O III] corresponds to the trough at ^ —1000 km s~^. 




Fig. 2. The H^-[Oiii] region of AKARI J1757+5907. (a) 
Continuum subtracted spectrum. The dotted line is the 
Fell template produced from I Zw 1. The Fell tem- 
plate is broadened and scaled at 5170 A. Note that the 
Fell template is clearly different from Fell emission line of 
AKARI J1757+5907 at 4600-4700 A and 5260-5400 A. (b) Fit 
of the spectrum. H/3, [Olll] doublet, and Fell AA4925,5020 
are fitted with a three Gaussians (cyan, blue and red lines), 
two sets of two Gaussians (green lines), and two Lorentzians 
(magenta lines), respectively. 
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Fig. 4. Top; Absorption troughs from He I A3890, 3190, and 
2946 as red, green, and blue histograms, respectively. Vertical 
bars represent the statistical uncertainties in the residual in- 
tensity. Bottom; Velocity-dependent column density deter- 
minations for He I* from the ADD (red histograms), C{v) 
(blue histograms), and power-law (green histograms) solu- 
tions. The velocity-integrated He I* column density values 
for the three methods and their associated statistical uncer- 
tainties are also listed. 



Fig. 5. Similar to Figure 4. Top; This shows the trough 
profiles for the Mgll doublet lines AA2796, 2804 (blue and 
red histograms, respectively) and their associated statistical 
uncertainties. Velocity-dependent column density determina- 
tions from the ADD (red histograms), C{v) (blue histograms), 
and power-law (green histograms) solutions. The velocity-in- 
tegrated column density values for the three methods and 
their associated statistical uncertainties are also listed. 
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Fig. 7. Reddening of AKARIJ1757+5907. Ordinate is a flux 
density corrected for the Milky Way extinction in units of erg 
cm~^ and abscissa is a rest wavelength. The rest 

optical spectrum is shown with the GALEX photometry (the 
most left cross) and continuum measurements (four crosses) 
from our spectrum. The power-law continuum reddened by 
the SMC-type extinction law with E{B — V) = 0.18 is shown 
with a dot-dashed lino. 



Fig. 6. Similar to Figures 4 and 5. Top: The trough pro- 
files for the Feu A2600 and 2587 lines as red and blue his- 
tograms, respectively. Their associated statistical uncertain- 
ties are shown as vertical slashes. Bottom: The AOD, C{v), 
and power-law column densities plotted as a function of ve- 
locity and the respective integrated values in red, blue, and 
green, respectively. 
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Table 1. Photometry of AKARI J1757+5907. 



Source 


Band 


nux density /magnitude 


GALEX 


NUV 


230.55 ± 8.03juJy 




9 


T R 1 m o re 
i_(J.i llldjg. 


POSS-I 


r 


15.0 mag. 


POSS-II 


9 


15.2 mag. 


POSS-II 


r 


15.3 mag. 


POSS-II 


i 


15.1 mag. 


2MASS 


J 


13.728 ±0.021 mag. 


2MASS 


H 


13.074 ±0.024 mag. 


2MASS 


Ks 


12.344 ±0.023 mag. 



*This value corresponds to (1.345 ± 0.047) X 10 '-^ erg s ^ cm ^ 
at 2267 A. 



Table 2. Properties of emission lines of AKAR1J1757+5907. 





Aobs 




FWHMtruc 


Line 


A 


z 


km s^^ 


H/3 


7872.0 


0.6189 


5080 


[Olll] 5008.24 


8088.3 


0.6150 


710 


[Olll] 5008.24 


8063.3 


0.6100 


480 



Table 3. Measured and model predicted column densities 



Measurement Method / Model parameter Nue \** Np^ ii(o)* ^Mg ii* 

Measurements 

Apparent Optical Depth 156.9 ±0.8 24.9 ±0.6 160.9 ±0.7 



Covering factor 
Power Law 








179.3 ±8.0 
180.2 ±3.2 


29.1 ±3.0 
30.4 ±3.7 


226.7 ±6.7 

602 ±246 








Models 








SED Z/Zq 


log Nh 




logn^ 








MF87 1.00 


20.82 


-2.15 


3.8 


176.5 


30.9 


1341.5 


Soft 1.00 


20.81 


-2.15 


3.8 


180.5 


32.4 


1364.6 


MF87 4.23 


20.66 


-2.17 


3.8 


172.3 


31.2 


1815.5 



*in units of 10^^ cm ^ 



Table 4. Properties of Measured Outflows to Date 



Object 



log LboI 

(ergs s~^) 



R 
(kpc) 



log Nh 
(cm^^) 



log Uh 



log Ek 
(ergs s~^) 



Reference^ 



AKARI J1757±5907 47.57 

QSO 1044+3656 46.84 

QSO 2359-1241 47.67 

SDSS J0838±2955 47.53 

SDSS J0318-0600 47.69 



>3.7 
1.7 ± 0.4 



3.2 tli 

-1.5 
-1.0 

5.9 ± 0.4 



3.3 



> 20.82 
20.84 ± 0.10 
20.56 ± 0.15 
20.80 ± 0.28 
19.90 ± 0.17 



> -2.15 
-2.19 ± 0.10 
-2.40 ± 0.15 
-1.95 ± 0.21 
-3.08 ± 0.05 



> 43.30 



44.81 



+0.09 



-0.11 

43.36 ± 0.27 
45.35 1q 22 
44.55 tVil 



> 70 
120 ± 25 
90 
300 
60 ± 20 



+35 
-20 
+210 
-120 



'Calculated using equation (3) assuming S2 = 0.2. The values for the last three objects are half of those founds in the reference due to the 
use of an improved estimate for Mq and E^. given by equation (3), over those given by equations (9) and (10) in Dunn et al. (2010). 
t 1-This Work, 2-Arav et al. (2010), 3-Korista et al. (2008) & Bautista et al. (2010), 4-Moe et al. (2009), 5-Dunn et al. (2010) 



